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A set of reduced Hall magnetohydrodynamic (MHD) equations are used to evaluate the stability 
of large aspect ratio current sheets to the formation of plasmoids (secondary islands). Reconnection 
is driven by resistivity in this analysis, which occurs at the resistive skin depth dr^ = S^'^^'^ y^LvA/'J, 
where Sl is the Lundquist number, L the length of the current sheet, va the Alfven speed, and 7 the 
growth rate. Modifications to a recent resistive MHD analysis [N. F. Loureiro, A. A. Schekochihin, 
and S. C. Cowley, Phys. Plasmas 14, 100703 (2007)] arise when collisions are sufficiently weak that 
drj is shorter than the ion skin depth di = c/upi. Secondary islands grow faster in this Hall MHD 
regime: the maximum growth rate scales as (di/L)^^-^^ S^^^'^^va/L and the number of plasmoids as 
{di/Ly^^^Sl^^'^^, compared to S]I^va/L and 5'^/^, respectively, in resistive MHD. 

PACS numbers: 52.35.Vd,52.22.Tn,94.30.cp,96.60.Iv 



I. INTRODUCTION 

The primary focus of modern magnetic reconnection 
research has been to explain why observed reconnection 
rates in high Lundquist number plasmas are much faster 
than predicted by the seminal Sweet-Parker theory!^ 
Important events such as the eruption of flares in the 
solar corona {Sl ^ 10^^) and the sawtooth collapse in 
magnetic fusion experiments {Sl ^ 10^) occur much 
faster than the resistive diffusion timescale predicted by 
the classical models.^ ^ Here Sl = ^t^Lva/ {(?Vi) is the 
Lundquist number based on the system size L, va is 
the Alfven speed, and r] the resistivity. The discrep- 
ancy between observed and predicted timescales has led 
to the point of view that fast reconnection cannot occur 
within the framework of resistive magnetohydrodynamics 
(MHD), which the Sweet-Parker theory is based on. 

Recent work has uncovered a fundamental flaw in the 
Sweet-Parker model when applied to high-Lundquist- 
number plasmas .'^'^II These studies have shown that the 
plasmoid instability, i.e., secondary island instability, 
which is excited when Sl ^ 10^, can signiflcantly en- 
hance the reconnection rate. In plasmoid dominated re- 
connection, the primary mechanism of reconnecting fleld 
lines is the formation of magnetic islands (flux ropes in 
These super-Alfvenic instabilities quickly grow, 
nonlinearly saturate, or coalesce, and convect out of the 
current sheet, carrying magnetic flux with them. Co- 
pious ejection of flux ropes have been observed at re- 
connection sites in solar flares^ as well as the Earth's 
magnetopause^^ and magnetotaiP^. 

Although secondary islands have been studied for a 
long time,^^ only recently have the scaling properties of 
the most unstable mode been established for a Sweet- 
Parker current sheet The primary insight has been to 
account for the Lundquist number scaling of the current 
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sheet width: (5sp = LSj^ ' . In conventional tearing 
mode theory, the current sheet width is taken to be con- 
stant, and the subsequent growth rate scales as Sl to 



a negative exponent: S^ ' or , for the const ant- 

ip and nonconstant-V^ regimes, respectively.^^ Accounting 
for the crucial feature that a Sweet-Parker layer becomes 
increasingly singular at high Sl^ the classical dispersion 
relation for tearing modes shows that the growth rate of 
the most unstable plasmoid scales as 5)^^ ^.1^ Thinning of 
the current sheet at high Sl plays a critical role, lead- 
ing to the surprising result that the plasmoid instability 
becomes increasingly unstable at higher Sl- 

These resistive MHD studies have provided a convinc- 
ing argument for the importance of plasmoids at high Sl^ 
but resistive MHD is not valid when is too high. Two 
common scenarios by which the resistive MHD approxi- 
mation breaks down are: (1) the current sheet width be- 
comes shorter than the ion skin depth di = c/ujpi^ (2) the 
resistive skin depth for the plasmoid dr, = Ssp \/ {va/ L)/"^ 
becomes shorter than the ion skin depth. Here 7 is the 
plasmoid growth rate. The resistive skin depth is the 
length scale at which plasmoids form and, thus, magnetic 
fleld lines reconnect. Scenario (1) is what is typically con- 
sidered Hall reconnection. If the current-sheet width is 
smaller than d^, a large aspect ratio current sheet (with 
Y-points) is no longer a viable equilibrium; it is typically 
(but not always, see Ref. TT) replaced by an X-point 
geometry, with an additional out of plane quadrupole 
fleld. Reconnection proceeds much faster in the X-point 
conflguration.^^ Starting from a Sweet-Parker current 
sheet, plasmoids can cause a cascade to the di scale 
because the current layer between plasmoids scales as 
3 Ssp/VN where N is the number of plasmoids. 

In this work, we consider scenario (2). Since the plas- 
moid instability is super-Alfvenic, 7 ^ va/L^ the resis- 
tive skin depth is much shorter than the current sheet 
width: drj/Ssp = (va/L)/^ <C 1. Thus, a Sweet- 
Parker current sheet that is thinning because of plas- 
moid production, or increasing Sl^ will always enter a 
regime where Hall effects alter the plasmoid instabil- 
ity [scenario (2)] before the current sheet becomes thin 
enough to alter the equilibrium [scenario (1)]. This in- 
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FIG. 1. Diagram of different reconnection scenarios in a 6i- 
Sl space. Also shown are the parameter range for several 
plasmas of interest: the magnetic reconnection experiment 
(MRX), Earth's magnetotail, the base of the solar chromo- 
sphere, and the solar corona. 



carried out by asymptotically matching solutions in the 
three layers (i) x <C (ii, (ii) dj^ <^ x <^ 6sp, and (iii) 
X ^ di. It is shown that the basic equations describ- 
ing the resistivity-driven Hall reconnection are similar to 
a collisionless reconnection problem solved by Mirnov, 
Hegna, and Prager in Ref. [2TJ and Fitzpatrick and Por- 
celli in Ref. 22, The physics is different because mag- 
netic field lines reconnect due to resistivity at the dr^ 
scale here, instead of electron inertia at the d^ scale in 
the collisionless problem. Nevertheless, the same bound- 
ary layer analysis can be used here. The properties of 
a Hall-plasmoid reconnection scenario are discussed in 
Sec. |V[ The analytic results are compared with linear 
and nonlinear simulations in Sec. IVll Section |Vl] also 
shows numerical solutions for the eigenfunctions, includ- 
ing the quadrupole out-of-plane magnetic field that ar- 
rises in the Hall-plasmoid regime. Results are summa- 
rized in Sec. I VIII 



termediate region, where (i^ <C (ii <C ^sp, provides a 
transition between plasmoid dominated reconnection in 
resistive MHD, where di <^ drj <^ 5sp, and conventional 
Hall reconnection, where 5sp <^ di. It is important to 
check that plasmoids continue to be unstable in this in- 
termediate regime, otherwise it could prevent the cascade 
to shorter scales. However, we find that this is not a con- 
cern because plasmoids are formed even more copiously 
in the Hall-plasmoid regime: the maximum growth rate 
scales as ^Yi/^o — ^\ 
moids as A^h — 



6/13^7/13 
q11/26 



and the number of plas- 
compared to "yYi/^o — S\l^ 
and A/r ^ 5^/^, respectively, in resistive MnD.!^ Here 
5i = di/L and F^ = 2va/L. 

The Hall-plasmoid region we are interested in is defined 

— 5/8 

by the boundaries dr, ^ di ^ (^sp, which implies Sj^ < 



]^ /2 

di < Sj^ • ^his last expression, either 7r or 7h can be 
used in d^ since both give 7/F0 5^^^ at di 



d^ for the 



fastest growing mode. Figure [T] shows these boundaries 
in SiSL space. Plasmas to the left of the vertical line 
at Sl = 10^ are considered stable to plasmoids. This 
result is based on numerical simulations which find that 
the plasmoid growth rate is sub-Alfvenic for Sl ^ 10^, so 
any islands convect out of the reconnection layer before 
growing to significant amplitudes. A large aspect ratio 
current sheet (Y-points) is often not a viable equilibrium 
for plasmas above the di = ^sp line. Plasmas below 
the di = dr^ line are in the conventional resistive MHD 
regime. The Hall-plasmoid region we are concerned with 
is a narrow part of the Si-SL parameter space at small 
Sl, but expands to cover a larger range of 5i for high- 
Sl plasmas such as the solar corona. Even for lower Sl, 
this region is important because it is always traversed as 
plasmoids cause a cascade from the resistive MHD regime 
to the Hall regime at shorter scales. 

Sectio n [ll| describes the reduced Hall MHD equations, 
and Sec. |III| the equilibrium, that are used in the linear 



II. REDUCED HALL MHD EQUATIONS 



In terms of the normalized variables 

a va CL 
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the single fluid equation of motion and the generalized 
Ohm's law can be written 



(St + V • V) V = J X B - VP, (2) 



and 



E + V X B = S-^J + di{J X B - VP), (3) 

respectively. The hats on normalized variables in Eqs. ([2| 
and ([3| have been omitted for notational convenience. 
Normalized variables will be used throughout Sees, [ill 
pl|and [lV| Here va = Po/a/^ttp is the Alfven speed and 
ujpi = ^/Aire'^n/mi is the ion plasma frequency. We will 
also use the Maxwell equations V-B = 0, VxE = — 9tB, 
and V X B = J. Equations Q and ([3| assume that pres- 
sure can be represented by a scalar, and Eq. ([3| neglects 
electron inertial scale physics {de 0). We assume that 
S~^ de, SO reconnection is driven by resistivity. It 
has been shown recently that in collisionless {S~^ <^ de) 
and weakly collisional {S~^ ^ de) regimes, mo re gen eral 
tensor descriptions of pressure can be required) ^"^ * * 

We assume incompressible flow, V • V = 0, and intro- 
duce the stream function 6, 



tearing mode analysis in Sec. IV The stability analysis is 



V = V^ X i + V^z, 

and flux function 

B = V?A X i + B^z. 



(4) 



(5) 
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In terms of the stream and flux functions, Eqs. ([2| and 
(|3| can be written as the fohowing set of four reduced 
resist ive-Hah MHD equations!^ 



(6) 



(7) 



dt^ = S-'V'^ + [(/), ^] + d,[i^, B,] + Eo, (8) 




dtB, = S-^V^B, + [c^,B,] + +d,[VV,V^]. (9) 

Equations (|6| and ^ are obtained from the curl of 
Eq. Q in the perpendicular (to z) and parallel directions, 
respectively. Likewise, Eqs. ([8| and (|9| are obtained from 
the curl of Eq. ([3|. Here Vq and Eq are constants, and 
[fid] = (^/ ^9) ' ^ is the Poisson bracket. Similar 
reduced Hall MHD equations are also derivable from gy- 
rokinetics, if a strong guide fleld is present 



III. EQUILIBRIUM 

To compare with the resistive MHD analysis of 
Loureiro et al.,^ we want to use the same Sweet-Parker 
equilibrium conflguration. Although we are using Hall 
MHD equations here, it is still possible to construct an 
equilibrium that does not depend on the ion inertial 
scale terms because the current sheet thickness is much 
larger than di. We choose the same linear flow proflle 
as Louriero: Vxo = — LoX, Vyo = Toy inside the current 
sheet {—Xo < x < Xo), Vxo = —Vo^ Vyo = above the 
current sheet {x > Xo), and Vx = Vo^ Vyo = below the 
current sheet {x < —Xq) . Here = 2a /L in normalized 
units {To = 2va/L in dimensional units). The associated 
stream function proflle is 



-Toxy, \x\ < Xo 

^o^oyi ^ 
To^oy^ ^ ^ 



(10) 



Assuming Vzo and B^o are constant, and looking for 
a steady-state ID solution of the form ipo = V^ol^), the 
only non-trival reduced equation is Eq. Q, which is 



1 dB,, 



xBo, 



(11) 



ToS dx 

for I a; I < Xq- Taking Byo{x = 0) = as the boundary 
condition, and matching the solution of Eq. ( '. 
±1 for \x\ > ±Xo, yields Louriero's ID equil: 



Byo 



Q;exp(— a;^/(5gp)^erfi(a;/(5sp), \x\ < Xo 



1, 

-1, 



x>Xo (12) 



in which erfl is the imaginary error function. Here 
^SP = ^Sp/<^ is the normalized Sweet-Parker width. The 



FIG. 2. Schematic drawing of three important scales in the 
Hall-plasmoid analysis: d^^, di, and Ssp- The sheared equi- 
librium magnetic field, Bo, and flow, Vo, profiles are also 
shown. 



transition point Xo is determined from the maximum of 



Byo, which gives Xo = 0.92(5sp and a = 1.85. 
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The detailed functional form of Eq. (12) will not be 



important in the inner region of the linear tearing mode 
analysis to follow. The critical feature is that the width of 
Byo is Ssp. Thus, the current sheet becomes increasingly 



singular as the resistivity decreases. Equation ( 12 ) shows 



that ^sp is a convenient parameter to use to normalize 
length scales. We choose the length a, which has been 
arbitrary until now, to be 



^SP, 



(13) 



so that 6sp = 1. The schematic drawing of the sheared 
magnetic fleld and flow proflles is shown in Fig. [2j along 
with the three scales, d^, di, and dsp, that will be im- 
portant in the tearing mode analysis to follow. 



IV. LINEAR GROWTH RATE 
A. Linearized equations 

We linearize Eqs. (|6|-(|9| according to 



(14) 



with analogous deflnitions for ^, Vz, and B^. We as- 
sume that the instability growth rate is much faster than 
the timescale for equilibrium flows into or out of the 
layer: 7 dx(l)oidy(j)o ^ T^, i.e., it is super- Alfvenic. In 
this case, the linear shear flow only contributes through 

Applying these assumptions, the linearized version of 
Eqs. (j6])-(j9| in the reconnection layer reduce to: 



7V^i = ikByoB^i, 



(15) 



7(9^ - k^)4., = ik[ByM - k^) - BU^i, (16) 



7[1 - - = ikByo4>i - diikByoB.i, (17) 
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j[l-dl{dl-k^)]B,^=ikBy,V,i 

+ diik[Byo{dl 



(18) 



In resistive MHD, tearing mode behavior is determined 
from just two equations: those for the the stream and flux 
functions [(16) and (17)]. In Hah MHD, four equations 
must be solved because the out of plane velocity and 
magnetic field are important at spatial scales of order di. 



B. Inner region 

In the inner region, x <C 1 (x <C ^sp in unnormalized 
units), so d'^ ^ k'^ and Byo — ax. Here, Eqs. (15)-(18) 
reduce to: 



^(1 — d^)ilJi = ix(j)i — idixBzi^ (19) 



^(1 -d^dl^ x^/g'^)Bzi = idixd'^tpi 



(20) 



Following Ref . E 
is 

TT 

in which 



the growth rate (in normalized units) 



dr^diG{g/di 



G{x) 



2 d,G{g/d,) g 
^/xT{l/4^x/4) 



2 r(3/4 + x/4)' 



(24) 



(25) 



and r(x) is the Gamma function. Since the growth rate 
analysis is based on g/di <C 1, we will use 



G{g/di 



i r(l/4) [J 
2 r(3/4) V di ■ 



(26) 



The general form of Eq. ([24|) was first obtained by 
Mirnov, Hegna, and Prager.^ Calculating the tearing 



mode growth rate from Eq. ( 24 ) requires the tearing sta- 
bility index. A', which is determined by tjji in the outer 
region. 



gdl(t)i= ixdli/ji 



(21) 



C. Outer region 



in which Vzi has been eliminated by putting Eq. ([15| into 
(18), and g = ^/{ak). 

Equations (Il9|-(l2T| are similar to Eqs. (41)-(43) of 



the collisionless tearing mode problem solved by Fitz- 
patrick and Porcelli.'^ In fact, making the substitutions 
de dr, and C/3 = V/3/(l + /3) ^ 1 in Ref. 122| yields 



identically Eqs. (19)-(|2T|. Fitzpatrick and Porcelli solve 



their equations using boundary layer theory by splitting 
the inner region into two parts: the narrowest part on 
the electron inert ial scale, x ^ d^ <^ di^ and a broader 
part on the ion inertial scale, <C (x ~ d^) <C 1. A 
growth rate solution is obtained by matching these two 
layer equations, as well as the ion inertial layer equation 
to the ideal MHD solution at large x. 



2 1 
A^x 



+ (9(x-2) 



where A' is the tearing stability inde:xP3 



A' 



MO) 



dt/ji 



dx 



0+ 



dipi 
dx 



(22) 



(23) 



Ref. 



The primary expansion parameter in 

g/{cf3di) < 1. 

Fitzpatrick and Porcelli's analysis can be carried over 
to solve the Hall MHD problem described by Eqs. ([l9|- 
(21 ) as long as the orderings dr^ <^di (in normalized 
units) are obeyed. In this case, the two parts that the 
inner region is split into are: the narrowest part on the 
scale of the resistive skin depth, x ^ drj <^ di^ and a 
broader part on the ion inertial scale, <C (x d^) <C 1. 
Here, the expansion parameter is g/di <^ 1, which will be 
checked a posteriori based on the growth rate solution. 



In the outer region x ^ di^ the out of plane magnetic 
and velocity perturbations decouple from the fiux and 
stream function equations because the di term in Eq. (17) 
is negligible. The term in Eq. (17) is also small, so 



the outer region is an ideal MHD solution. Here, Eq. (17) 
reduces to = —i^/{kByo)ipi^ and Eq. (16) to 



^1 



yo 



(27) 



Equation |27) describes the same outer region as the 
ysis of Loureiro et al,^ who solved it per- 

'x\ < Xo 



resistive ana 

turbatively by matching solutions from the 



{dl > k^) and \x\ 
The solution in the 



> Xo {d^. <C k'^) asymptotic limits. 
\x\ < Xo limit is 



^t{x) = ctByo{x)^a 



~-Byo{x) I 



dz 



(28) 



Taking Byo to be continuous at x = and using Byo — 
ax for X <C 1 yields = — a?/^i(0). The solution of 



Eq. (27) in the \x\ > Xo region is ip^ = exp(^/cx). 



Matching the first derivative of these two solutions ai x = 
Xo, using Byo{^Xo) = ±1 and By^{±Xo) = yields = 
<^'0i(O) exp(/cXo)//c. Matching the solutions themselves 
yields = ±ailji{0)/k. Thus, 



±- 



^V^i(Q) r Byo[l T k jl^^ dzB-^] \x\ < Xo 



k \ exp[/c(xo ^ x)\ 
which yields the tearing stability inde:jP 



\X\ > Xo 



'(29) 



A' 



2a2 
1^' 



(30) 
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FIG. 3. Growth rate calculated from Eq. (|3l| using fixed 
ion skin depth 6i = 10~^, and three values of the Lundquist 
number Sl = 10^ 10^ and 10^ 



6i^l X 10-^ 
(5, = 3 X 10"^ 
1 X 10"^ 



-1/2 




FIG. 4. Growth rate calculated from Eq. (31) using fixed 



Lundquist number Sl — 10 , and three values of the ion skin 
depth 5^ = 1 X 10"^3 X 10"^ and 1 x 10"^. 




FIG. 5. Contour plot of the peak growth rate in the 
resistive-plasmoid {'^ /To — S]^"^) and Hall-plasmoid (7/ro — 

3:6/13 o7/13n 

d/ bj^ ) regions. 



Eq. (31) is negligible. Here, the growth rate is given 

by 



J_ _ (aCg)^/^ 4/5 c2/5ol/5 



(33) 



For large k, the first term in Eq. (31) is negligible. Here, 
the growth rate is given by 



1. 
To 



(34) 



The Kf^^^ and n scalings predicted by Eqs. (33) and 
(34) match well with the numerical solutions of (31). 



Figures |3] and |4] show that a broad range of wavenum- 
bers can be unstable, but the most unstable will dom- 
inate the reconnect ion process. The wavenumber of 
the most unstable mode can be obtained by equating 
Eqs. (33) and (34), which yields 



V. HALL PLASMOID RECONNECTION 



2^^^1/13^11/26 



(35) 



In terms of unnormalized, i.e., dimensional, variables, 
the dispersion relation from Eqs. (24), (26) and (30) is 



\ 5/2 5/2rl/2 



(8a)i/25y4r, 2v^ 



in which n = kL^ Si = di/L^ and 



' Sl^^ = 0, (31) 



1 r(i/4) 

2 r(3/4) 



1.48. 



(32) 



Figure [3] shows j/To from a numerical solution of Eq. (31 ) 
for fixed Si = 1 x 10~^, and three values of Lundquist 
number, Sl = 10^, 10^ and 10^. Figure |4] shows the 
solution for fixed Sl = 10^, and three values of the ion 
skin depth Si = 1 x 10"^3 x 10"^ and 1 x 10""^. 

Figures [3] and 4] show that 7/ro has a power-law de- 
pendence with different slopes in the large and small 
wavenumber limits. For small a^, the second term in 



Putting Eq. (|35| into either Eq. ([33|) or Eq. (|34|) provides 
the growth rate of the most unstable plasmoid 



7max/ro 



3 2^6/13^7/13^ 



(36) 



Since the number of plasmoids in the chain may be esti- 
mated to be :^ /^max/27r, Eq. (35) shows that the num- 
ber of plasmoids formed for a given Si scales slightly more 
rapidly with Lundquist number in the Hall-MHD regime 
than the resistive-MHD regime (where A^max — ^i!^)- 
The growth rate of the most unstable mode also grows 
more rapidly with Sl (7max/ro — S]^^^ in the resistive 



case). Equation (35) and Fig. [4] show that A^max is quite 
insensitive to the ion inertial length, being proportional 
to Sy^^. However, Eq. (35) and Fig. [3] show that A^max 
is sensitive to Lundquist number. In contrast, the max- 
imum growth rate is sensitive to both Si and Sl- Fig- 
ure [5] shows the peak growth rate as a function of Si and 
Sl in the resistive-plasmoid and Hall-plasmoid unstable 
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FIG. 6. Lin ear growth rate calculated numerically from 
Eqs. (15)-(18) using a fixed Lundquist number Sl — 10^, and 
three values of the ion skin depth = 1 x 10"^ 3 x 10"^ 
and 1 X 10"^. 



regions. The linear stability of plasmoids in the conven- 
tional Hall regime, di > Ssp^ has yet to be worked out. 

Finally, we check the assumptions made during our 
analysis: drj <^ di <^ Ssp and g/di <C 1. The drj <^ di <^ 



-5/8 



< S,^ < S] 



-1/2 



For Si 



Ssp assumption requires Sj^^ 

10~^, this implies 10^ ^ ^ 10^, which determines the 
Sl limits in Fig. [sj Likewise, for Sl = 10^, this requires 
10~^ ^ Si ^ 10~^, which determines the Si limits in 
Fig. |4j The g/di {in normalized units) criterion can 
be written 



a K,5i r„ 



< 1. 



(37) 



-5/8 



For the most unstable mode, this reduces to Si> S\^ 
which is the same requirement already obtained from the 



dj^ ^ di condition. Equation (|37|) is more restrictive for 



1^ 7^ /^max, but wc are primarily interested only in the 
most unstable mode here. Equation (30) assumed n <C 



5' 



1/2 



which sets the maximum n that can be considered 



in Figs. [3] and [4j but does not affect the most unstable 
mode (this will be discussed more in Sec. VIA). We also 



assumed that ^ /Tq ^ 1, which sets the minimum n that 
can be considered in Figs. |3] and [4j 



VI. NUMERICAL RESULTS 



A. Numerical Linear Eigenmode Solutions 



The analysis of Sec. [IV|is based on the ordering <C 
di <^ ^SP7 but Fig. [5] shows that the range of di satis- 
fying this can be narrow, especially at lower Sl- Thus, 
it is useful to check numerically the linear growth rate 
formulas obtained in Sec. [Vl To do so, we solve numer- 
ically the four linearized reduced Hall MHD Eqs. (15)- 
(18). Figure [g] shows the plasmoid growth rate using a 
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FIG. 7. Eigenfunctions for the perturbed fields V^i, 01, Bzi^ 



and Vzl calculated numerically from Eqs. (15)-(18) using 
Lundquist number Sl — 10^, and two values of the ion skin 
depth: (a) a resistive MHD case with = 0, and (b) a Hall 
MHD case with = 3 x 10~^. All eigenfunctions are normal- 
ized by their maximum amplitude. 



fixed Lundquist number S'^, = 1 x 10^, and three val- 
ues of the ion skin depth Si = 1 x 10~^, 3 x 10~^, and 
1 X 10~^. These curves can be compared to the analytic 
results from Fig. [4] 

For Sl = 1 x 10^, the boundary between the resistive 
plasmoid and Hall-plasmoid regimes {dr^ = di) occurs at 

Si ^ S~^^^ = 1 X 10-^ 



at Si 



S] 



-1/2 



The boundary di = ^sp occurs 
1 X 10~^. Figure 6| shows that for 



di = 1 X 10 ^ — dn^ the resistive MHD"^ scaling with 

is found. As Si increases through the Hall plasmoid 
regime, this steepens, in agreement with the theory. The 
Hall MHD scaling, from Eq. (33), is found for the 

larger Si. 

The numerical results of Fig. |6] show that stabilization 
occurs quickly after the maximum k considered in the an- 
alytic theory is exceeded. The A' solution used in the an- 
alytic analysis, Eq. (30), assumed kS^p <C 1 (/^ <C Sj^ ). 
A more complete A^olution that can account for this 
stabilization is A' ^ 2a^S]l'^/i<i + 0.57 - lAln^S^^ (see 
Ref . HI) . Here we are primarily interested in the growth 
rate of the most unstable mode, which is accurately pre- 
dicted by the approximate expression for A^ 

Figure [7| shows profiles of the eigenfunctions ^/^i , ^1 , 
Bzi^ and Vzi for Sl = 10^ and two values of Si. Each 
eigenfunction is normalized by its maximum amplitude. 
Panel (a) shows a resistive MHD solution, where the 
ion skin depth was taken to vanish Si — 0. Here, the 
out of plane magnetic field and flow are zero, and the 
stream and flux functions of conventional MHD tear- 
ing modes^ are found. Panel (b) shows a case with 
^ • = 3 X 10~^, which, according to Fig. [sj is just into the 
Hall MHD regime (d^ < < ^sp) for this Sl^ Figure [t] 
shows that the stream and flux functions have similar 
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ky /2t: 

FIG. 8. Contour plot of constant Bzi surfaces for a single 
period of the plasmoid chain. Bzi has a quadrupole form, 
which is characteristic of the out-of-plane magnetic field in 
Hall reconnection. Here Sl — I ^ 10^, and = 3 x 10~^ 



profiles in both the Hall and resistive MHD regimes, but 
that the out of plane magnetic field and flow velocity be- 
come important for spatial scales shorter than di in Hall 
MHD. 

Figure [8] shows a contour plot of constant Bzi sur- 
faces over a single period of the plasmoid chain. The 
out-of-plane magnetic field generated by the presence 
of magnetic islands for this Hall MHD regime has a 
quadrupole form. Similar quadrupole fields have been ob- 
served in Hall MHD simulations of steady reconnection.^ 
However, these are typically seen when di > Ssp^ and 
the equilibrium has an X-point configuration. Here 
quadrupole fields are formed due to the presence of plas- 
moids in an extended current-sheet equilibrium. 



B. Evidence of the Instability in Nonlinear 
Simulations 

The prediction that Hall effects modify the plasmoid 
instability when < di can also be tested using more so- 
phisticated codes that solve the full nonlinear Hall MHD 
equations. Here we employ the same simulation setup 
of two coalescing magnetic islands as in Refs. "9^ and TP, 
Length scales are normalized to the length of the simu- 
lation box (L = 1), which is square. Initially, the cur- 
rent sheet between the two islands is thicker than the 
Sweet-Parker width. As the simulation proceeds, the cur- 
rent sheet gradually becomes thinner until it reaches the 
Sweet-Parker width. To obtain the linear growth rate 
of plasmoids, we integrate B^ at the central part of the 
current sheet along x = 0, from y = —1/4 to 1/4 at each 

time step. The magnitude of f{t) = J^y^Bl{t)dy re- 
mains small before onset of the plasmoid instability, and 
increases abruptly after the onset. We fit ln(/) during 
the period that the plasmoid grows in accordance with 



Si 

FIG. 9. Linear growth rate of plasmoids obtained from a 
numerical simulation of a full nonlinear set of Hall MHD 
equations. Here Sl = x 10^, for which di — dr^ when 
= 1.8 X 10~^, and di = fep when = 1 x 10~^. 



ln(/) 2^t + c. The 7 so obtained is a good measure of 
the linear growth rate of the fastest growing mode. 

Figure [9] shows the linear growth rate inferred from 
nonlinear Hall MHD simulations using /Sl = 1 x 10^, and 
5i between 1 x 10"^ and 1 x 10"^. For 5^ = 1 x 10^, 
the Hall MHD regime is expected when di > d^, which 

implies 5i > Sl^^^ = 1.8 x IQ-"^. Figure [9] shows that the 
transition between the resistive-plasmoid regime, where 
the growth rate is independent of (5^, and the Hall- 
plasmoid regime, where the growth rate depends on Si, 
occurs around Si 2 x 10~^. This agrees well with the 
analytic expectation. The analytic theory is not expected 
to hold for Si > S^^^'^ = 1 x 10~^, due to changes in the 

equilibrium. It is difficult to confirm the predicted S^^^^ 
scaling in the Hall-plasmoid regime because of uncertain- 
ties associated with inferring the linear growth rate from 
nonlinear simulations. However, confirming that Hall ef- 
fects are important when the ion skin depth is shorter 
than (i^, rather than Ssp^ provides an important consis- 
tency check with the analytic theory. The Hall-plasmoid 
regime is expected to include a much broader range of Si 
for high-^L plasmas, such as the solar corona, but, un- 
fortunately, computational resources limit the maximum 
Sl obtainable in nonlinear simulations. 



VII. SUMMARY 

There are two scale lengths that, when short enough, 
can cause the conventional resistive-MHD tearing mode 
theory to break down: the current sheet width (^sp)^ and 
the resistive skin depth {dr^). When Ssp ^ di, the equilib- 
rium magnetic configuration is expected to change. An 
X-point geometry is typically seen in this regime, which 
is often referred to as the Hall reconnection regime. How- 
ever, if the plasmoid instability is present, the resistive 



8 



skin depth is much shorter than the current sheet width: 
dr] ^ ^^sp- Thus, Hah effects modify the plasmoid in- 
stabihty properties at lower Sl than is required to cause 
a change in the equihbrium configuration. As plasmoids 
cause a cascade to shorter scales, this means that the 
Hall-plasmoid regime is always realized before the scales 
become short enough to modify the equilibrium. Out of 
plane velocity and magnetic field components arise in the 
Hall-plasmoid regime, with the magnetic field having a 
quadrupole form. 

Plasmoids grow more rapidly with Sl in the Hall- 
plasmoid regime than the resist ive-plasmoid regime. 
The growth rate of the most unstable mode scales as 
d^^^^S]^^^^ and the number of plasmoids as S^^^ S]^^"^^ . 
The corresponding resistive-MHD scalings are and 

S^^^. Thus, the cascade to shorter scales that the plas- 
moid instability causes is expected to continue through 
the Hall-plasmoid regime until the relevant current-sheet 



width (between plasmoids) becomes shorter than di. At 
this point, the Sweet-Parker equilibrium, characterized 
by an extended current-sheet, will be superseded by a 
Hall-MHD equilibrium. 



ACKNOWLEDGMENTS 

The authors gratefully acknowledge conversations with 
Dr. Brian Sullivan and Dr. Will Fox. This research was 
supported in part by an appointment to the U.S. Depart- 
ment of Energy Fusion Energy Postdoctoral Research 
Program administered by the Oak Ridge Institute for 
Science and Education (S.D.B.), and DOE Grant No. 
DE-FG02-07ER46372, NSF Grant Nos. ATM-0802727, 
ATM-0903915, and AGS-0962698, and NASA Grant Nos. 
NNX09AJ86G and NNX10AC04G. 



^ P. A. Sweet, Nuovo Cimento Suppl. 8, 188 (1958). 
2 E. N. Parker, Astrophys. J., Suppl. Ser. 8, 177 (1963). 
^ A. Bhattacharjee, Annu. Rev. Astro. Astrophys. 42, 365 
(2004). 

^ M. Yamada, R. Kulsrud, and H. Ji, Rev. Mod. Phys. 82, 
603 (2010). 

^ N. F. Loureiro, A. A. Schekochihin, and S. C. Cowley, 

Phys. Plasmas 14, 100703 (2007). 
^ R. Samtaney, N. F. Loureiro, D. A. Uzdensky, A. A. 

Schekochihin, and S. C. Cowley, Phys. Rev. Lett. 103, 

105004 (2009). 

^ A. Bhattacharjee, Y.-M. Huang, H. Yang, and B. Rogers, 

Phys. Plasmas 16, 112102 (2009). 
^ P. A. Cassak, M. A. Shay, and J. F. Drake, Phys. Plasmas 

16, 120702 (2009). 
^ Y.-M. Huang, and A. Bhattacharjee, Phys. Plasmas 17, 

062104 (2010). 

^° L. Nei, K. Germaschewski, Y.-M. Huang, B. P. Sullivan, 
H. Yang, and A. Bhattacharjee, Phys. Plasmas 17, 052109 
(2010). 

Y.-M. Huang, A. Bhattacharjee, and B. P. Sullivan, Phys. 
Plasmas 18, 072109 (2011). 

D. A. Uzdensky, N. F. Loureiro, and A. A. Schekochihin, 
Phys. Rev. Lett. 105, 235002 (2010). 

W. Daughton, V. Roytershteyn, H. Karimabadi, L. Yin, 

B. J. Albright, B. Bergen, and K. J. Bowers, Nat. Phys. 7, 
539 (2011). 

S. L. Savage, D. E. McKenzie, K. K. Reeves, T. G. Forbes, 
and D. W. Longcope, Ap. J. 722, 329 (2010); R. Liu, J. 
Lee, T. Wang, G. Stenborg, C. Liu, and H. Wang, Ap. J. 
Lett. 723, L28 (2010). 

C. T. Russell, and R. C. Elphic, Space Sci. Rev. 22, 681 
(1978). 

L.-J. Chen, A. Bhattacharjee, P. A. Puhl-Quinn, H. Yang, 
N. Bessho, S. Imada, S. Miihlbachler, P. W. Daly, B. Lefeb- 
vre, Y. Khotyaintsev, A. Vaivads, A. Fazakerley, and E. 



Georgescu, Nat. Phys. 4, 19 (2007). 

H. P. Furth, J. Killeen, and M. N. Rosenbluth, Phys. Fluids 
6, 459 (1963). 

B. Coppi, R. Galvao, R. Pellat, M. Rosenbluth, and P. 
Rutherford, Sov. J. Plasma Phys. 2, 533 (1976). 
L. S. Shepherd and P. A. Cassak, Phys. Rev. Lett. 105, 
015004 (2010). 

W. Daughton, V. Roytershteyn, B. J. Albright, H. 
Karimabadi, L. Yin, and K. J. Bowers, Phys. Rev. Lett. 
103, 065004 (2009). 

V. V. Mirnov, C. C. Hegna, and S. C. Prager, Phys. Plas- 
mas 11, 4468 (2004). 

R. Fitzpatrick and F. Porcelh, Phys. Plasmas 11, 4713 

(2004) . 

R. Fitzpatrick and F. Porcelh, Phys. Plasmas 14, 049902 
(2007). 

M. Hosseinpour, N. Bian, and G. Vekstein, Phys. Plasmas 
16, 012104 (2009). 

R. Fitzpatrick, Phys. Plasmas 17, 042101 (2010). 
2^ L J. D. Craig, and P. G Watson, Phys. Plasmas 12, 012306 

(2005) . 

2^ A. A. Schekochihin, S. C. Cowley, W. Dorland, G W. 
Hammett, G G. Howes, E. Quataert, and T. Tatsuno, As- 
trophys. J. Suppl. Ser. 182, 310 (2009). 
The constant a is 1.31 in Ref. |5] because it uses 6cs — 
\/c^T\l (47rro) to characterize the current sheet thickness. 
Here we use fep, which is v^fcs- 

Parts of Ref.|22|have been debated recent ly.^*^ However, 
this debate is concerned with determining if the electron 
pressure model used in|22]is appropriate for modeling col- 
lisionless tearing modes. This concern does not affect the 
Hall MHD problem we are interested in here. 
T. Terasawa, Geophys. Res. Lett. 10, 475 (1983); B. N. 
Rogers, R. E. Denton, J. F. Drake, and M. A. Shay, Phys. 
Rev. Lett. 87, 195004 (2001). 



